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We investigate the effects of mobility and density on the amplitude death of coupled oscillators 
in metapopulation networks, wherein each node represents a subpopulation with any number of 
mobile individuals. We perform stochastic simulations of the dynamical reaction-diffusion processes 
associated with the Landau-Stuart oscillators in scale-free networks. Interestingly, we find that, with 
increasing the mobility rate or density, the system may undergo phase transitions from incoherent 
state to amplitude death, and then to frequency synchronization. Especially, there exists an extent 
of intermediate mobility rate and density leading to global oscillator death. In addition, we show this 
nontrivial phenomenon is robust to different network topologies. Our findings may invoke further 
efforts and attentions to explore the underlying mechanism of collective behaviors in metapopulation 
coupled systems. 

PACS numbers: 89.75.-k, 05.45.-a 
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I. INTRODUCTION 

The complex dynamics of coupled nonlinear oscillators 
has gained significant research interest in many branches 
of science and technology jl| , ranging from the modeling 
of biological rhythms in the heart [2|-y , nervous system 
@, HI: intestine [7j and pancreas @. One of the most in- 
triguing emergent phenomena is the so-called amplitude 
death (AD), a phenomenon that coupled oscillators stop 
oscillation and settle to a quenched steady state. Since 
the initial work by Bar-Eli [9|, AD has been observed in 
various real systems, such as the Belousov-Zhabotinsky 
reaction system relativistic magnetrons [ll|, and 

synthetic genetic networks p"3 - [l5j . to list just a few. 
The physical mechanism underlying AD has also been 
extensively investigated [iGl - fTgj . For coupled oscillators 
on a regular lattice, mechanisms leading to AD may in- 
clude oscillator mismatch p, IgI. delayed coupling [20l - 
[22|, asymmetric coupling p5-l27|, ti me-dependent cou- 
pling [28], conjugate coupling [23,|30[, nonlinear (nondif- 
fusiye) coupling [3l| , linear augmentation [sll , and so on 
[33l - [35| . See a recent review [l| for more details, which 
gives the deep insight and many clues for understanding 
AD induced by different scenarios. 

On the other hand, AD in complex networks of coupled 
oscillators has also attracted much attention [31, 36 43]. 
It has been shown that network topologies, being reg- 
ular (36l-[40j. random 41], small- world [i^] or scale- free 



[431 ] ■ may influence AD in a nontrivial way, by enhanc- 
ing or suppressing it with varying network randomness 
or heterogeneity. Nevertheless, previous studies on AD 
in complex networks only deal with the case of immo- 



bile individuals and each network node is occupied by 
one single individual. Very recently, the metapopulation 
network model which incorporates local node pop- 
ulation, mobility over the nodes, and a complex network 
structure, has drawn intensive attention. This model has 
been successfully exploited in different contexts, includ- 
ing epidemic spreading [45'~'48'] , biological pattern forma- 
tion , 49 , 50], chemical reactions 51], population evolu- 
tion Im], and many other spatially distributed systems 
[53I |54| . It is shown that the density and the mobility of 
the individuals could have drastic impact on the emer- 
gence of collective behaviors in general [3, [E^ , and par- 
ticularly, density and mobility induced synchronization 
of coupled oscillators has been reported [H, [13] ■ There- 
fore, one may ask: How would the mobility and the den- 
sity would influence AD in metapopulation networks of 
coupled oscillators? 

In the present work, we study the emergence of AD 
in a metapopulation model which incorporates mobility 
over a complex network together with local interactions 
of the individuals at the network nodes. We find that, for 
small mobility rate and lower density, the metapopula- 
tion displays incoherent state, for intermediate mobility 
rate or density, all the oscillators spontaneously cease 
their oscillations, but for large ones, the death state can 
be eliminated and synchronize to a common frequency, 
overcoming the disorder in their natural frequencies. Fur- 
thermore, we show this nontrivial phenomenon is robust 
to different network topologies. 



II. 



MODEL DESCRIPTION 
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We consider a system of N distinct subpopulations la- 
beled jjL, each corresponding to a network node. The den- 
sity p of the metapopulation is given by p = T^^=i 



2 



where iV^ is the number of individuals in node /z. Individ- 
uals inside each subpopulation run stochastically through 
the paradigmatic Landau-Stuart oscillators whose 
dynamics is described by: 

z^{t) = z^it){r-\z^\' + ^u;^) + Ki{z)^-z^it)){l) 

where z^{t) is the complex amplitude of the (j — 
l,...,Nf_,) oscillator in the /i-th node, r > denotes a 
homogeneous oscillation growth rate, K is the coupling 
strength, and 

side node fi. ujj denotes the natural frequency of the 
jth oscillator within node fi, which is picked up from a 
certain distribution. In the present work, we adopt a lin- 
ear distribution — 7 < w < 7 with 7 a constant. Without 
coupling, the trajectory of each single oscillator will settle 
to a limit cycle with frequency and radius = r. 

With sufficiently strong K and wide distribution of w^, 

and 



z)^ = ^fj^i Zi{t) is the mean field in- 



an 



the oscillators pull each other off their limit cycles, 
collapse into a state of zero amplitude = 0, i.e., 
AD state. 

The above equation actually defines the "reaction" 
process that governs the temporal behavior of each in- 
dividual inside the metapopulation nodes. We now as- 
sume that the individuals can diffusion randomly among 
the nodes. The system evolves in time according to the 
following rules 4J]. We introduce a discrete time step 
T representing the fixed time scale of the process. The 
reaction and diffusion rates are therefore converted into 
probabilities. In the reaction step, all the individuals are 
updated in parallel according to EqUJ After that, dif- 
fusions take place by allowing each individual to move 
into a randomly chosen neighboring node with probabil- 
ity Dt, where D denotes the mobility rate. If not other- 
wise specified, the parameters are N = 1000, r = 0.001, 
r = 0.4, 7 = 3.0, and K = 10. We choose the mobil- 
ity rate D and the density p as main control parameters. 
Each simulation plot is obtained via averaging over 20 
independent runs. 



III. RESULTS AND DISCUSSION 

To begin, we consider scale-free networks generated by 
using the Barabasi-Albert (BA) model [s^ with power- 
law degree distribution p{k) ^ k~^. We fix p = 10 (thus 
we have totally Np — Np individual oscillators) and vary 
D to investigate how the oscillators evolve in time. Ini- 
tially, the oscillators are homogeneous distributed among 
the nodes. If diffusion is absent {D = 0), each node will 
stay in an incoherent state for the above-mentioned pa- 
rameters. If D is small, each node still remains incoher- 
ent, as shown in Fig. 1(a) for D = 0.1, where typical time 
series of Re{Zj{t)) for several randomly-chosen oscillators 
within a random chosen node are depicted. For moderate 
diffusion rate D = 1.0 as shown in Fig. 1(b), the time se- 
ries eventually collapse into I ~ 1 , which indicates the 



D, say D = 10 in Fig. 1(c), the oscillators are locked to a 
synchronized state. Therefore, we observe an interesting 
mobility induced transition from incoherence to AD and 
then to synchronization in our metapopulation oscillator 
network model. In addition, we have also perform sim- 
ulations with fixed D and varying p. Similar behaviors 
(time series not shown) are found, i,e, the system will 
also undergo a density induced transition from incoher- 
ence to AD and then to synchronization with increasing 
p given D is not too small. 
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FIG. 1: (Color online) Time series of the real parts of the 
complex vector z^{t) for several random oscillators at D = 0.1 
(upper panel), D — 1.0 (middle panel) and D = 10 (lower 
panel). All the networks have fixed average network degree 
(k) — 6 and N = 1000. Other parameters are p = 10 and 
K = 10. 

To quantitatively manifest the collective oscillation in- 
tensity of the system, we choose the normalized mean 
"incoherent" energy E and "coherent" energy W as the 
characteristic variables, which are defined as 



M' as 



E 



W = 



Npr 



(2) 



occurrence of global AD. However, for sufficiently large Here the brackets (•) denotes averaging over time and 
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[■] stands for averaging over 20 different network real- 
izations for each D and p. A large E means relatively 
large average oscillation amplitudes, while a larger W 
implicates more synchronous oscillations [42]. In the in- 
coherent state, W will be nearly zero but E can be large. 
For AD state, both E and W should be nearly zero. In 
the synchronized state, E and W both have noteworthy 
nonzero values and they should be equal for complete 
synchronization. 

The dependences of E and W on D for fixed p — 10.0 
and p for fixed D — 0.4 are depicted in Fig[21 wherein the 
curves can be divided into three stages. In stage 1, E de- 
creases monotonously with increasing D(p), approaching 
nearly zero at a certain value of D{p). In stage 2, both 
E and W remain nearly zero, demonstrating the occur- 
rence of global AD. Both E and W increase with D{p) 
in stage 3, corresponding to the synchronized state. In 
accordance with Fig.l, the transition from incoherence to 
global AD and then to synchronization is clearly demon- 
strated. It seems that increasing density plays a similar 
role to increasing the mobility. Of particular interest, 
there exists a notable region (stage 2) for both D and p 
where the system remains stably in the AD state. 
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FIG. 2: (Color online) The "incoherent" energy E and "coher- 
ent" energy W of the limit cycles of individual oscillators as a 
function of D in panel (a), and p in panel (b). Three stages 
are presented as indicated by 1, 2, and 3. 

To further get a global picture, we have perform ex- 
tensive simulations to obtain the phase diagram in the 
D p plane. This is shown in Fig. 3, where the con- 
tour plots of E and W are drawn. As expected, when 
both D and p are small, the system shows the incoherent 
state (the light region 1 of Fig. [31[a)). When D and p 
are increased crossing certain critical values, the system 
finally evolves into global AD state (the dark region 2 of 
Figl3l[a)). For sufficiently large D and p, however, the 
system may undergo a phase transition again from the 
death phase to synchronized one, as shown in the light 
region 3 of Fig. [2Ka). Thus as we keep one of the param- 
eters fixed {D or p ) and increase the other, we observe 
transitions from the incoherent (region 1) to the ampli- 
tude death (region 2) and then to frequency synchronized 
state (region 3). Note that the 1 — > 2 transition is not 
observable with increasing p if the mobility rate D is too 
small, and similarly, the 2^-3 transition cannot happen 
if the density p is too low. 

So far, the results are all for scale-free coupled net- 
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FIG. 3: (Color online) Panels (a) and (b) correspond to phase 
diagram E{D,p) and W{D,p) of the metapopulation model 
respectively, both on a synthesized 1000-node BA network 
with (fc) — 6. 



works. One may wondering whether the interesting find- 
ings above is sensitive to the network topology or not. 
Thus, we have also performed similar studies on other 
types of networks, e.g., the small world network and ran- 
dom network. The phase diagrams are shown in Fig|4] 
(a,b) for small-world network and (c,d) for random net- 
work, respectively. Apparently, the qualitative behaviors 
are the same as those observed in scale-free networks. 
The only difference is that the boundaries between dif- 
ferent phase regions are slightly shifted. 
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FIG. 4: (Color online) Panels (a) and (b) correspond to phase 
diagram E{D,p) and W{D,p) of the metapopulation model 
on 1000-node small-world networks respectively, (c) and (d) 
for random networks. Other parameters are the same as Fig 
3. 
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IV. CONCLUSION 

In summary, we have studied the collective dynamics 
of coupled limit cycle oscillators defined on metapopu- 
lation networks, where different subpopulations are con- 
nected by fluxes of individuals. By extensive numerical 
simulations, we show that mobility and density play non- 
trivial roles on the collective behavior of the system, by 
demonstrating an interesting type of density or mobility 
induced transition into and out from oscillation death 
behavior. On one hand, intermediate mobility rate and 
density can induce global oscillator death. On the other 
hand, large mobility rate and density can also eliminate 
oscillator death and lead to synchronization. In addition, 
we find that this nontrivial phenomenon is robust to the 
network topology. The underlying mechanisms for these 
interesting results are still open. An analytical theory 
would surely be greatly helpful, however, it is not avail- 
able at the current stage and should certainly deserve 



more study. Since many real-life networks (cellular net- 
works, protein networks, gene networks, etc.) inevitably 
involve variances in both the mobility and the density, 
and their collective dynamics could be modeled by cou- 
pled oscillators, these results may find a variety of ap- 
plications. Our study may also stimulate further investi- 
gation on the emergent coherence of metapopulations of 
coupled oscillators. 
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